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Q ■ Abstract 

It is discussed the possibility of describing the Doppler acceleration 
towards the Sun measured by the Pioneer space-craft above the he- 
liocentric distance of 20 AU, known as the Pioneer anomaly, by con- 
sidering an expanding background described by an expanding locally 
anisotropic (ELA) metric. This metric encodes both the standard lo- 
cal Schwarzschild gravitational effects and global universe expansion 
\^ • effects allowing simultaneously to fine-tune other gravitational effects 

. at intermediate scales, such as local dark matter and/or gravitational 

interactions corrections. Taking the Sun as central mass and the 
CSJ , heuristic ansatz for the dimensionless functional exponent parameter 

' a = 3 -|- a2(ri)(l — J7q)^ with a linear ansatz Q!2(ri) = (22^''^' -|- ttj^'"'^^ ri 

above 20 AU, being Uq the Schwarzschild gravitational potential of the 
Sun, the Pioneer effect is fully accounted for as due to the blue-shift of 
the background space-time described by this metric. For compatibility 
' with orbital motion within the Solar system it is further considered 

%-{ . a branch ansatz for which a2{ri) = below 20 AU. This construc- 

tion has the advantage of reducing the physical acceleration correction 
above 20 AU down to ~ 10~^'^ ■ms~^ outwards from the Sun, hence 4 
orders of magnitude below the measured Pioneer acceleration such that 
orbital corrections to Neptune and Pluto are less notorious than when 
the full Pioneer acceleration is considered to be a gravitational acceler- 
ation. It is also discussed the effective mass-energy density distribution 
above 20 AU, the equation of state (—1 < ujr,uje < —0.726) due to the 
background corrections and the respective finite non-null total mass is 
computed by considering an upper cut-off for the functional parameter 
a. Such effective extended configuration is tcniptatively interpreted as 
a cold dark matter extended distribution. 
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1 Introduction 



The Pioneer space-crafts have been launch from Earth on 1972 and 1973 
travelhng outwards of the Solar system |l|2j . In between the heliocentric 
distances of 20 AU and 70 AU, travelling at an approximately constant ve- 
locity of Vp = 12.2 X 10'^ ms~^, it was detected an unmodeled Doppler shift 
known as the Pioneer effect. If totally attributed to a physical acceleration, 
this effect corresponds to a radial acceleration (towards the Sun) of [5] 

Op = -(8.74 ±1.33) X 10"^° ms"2 . (1) 

Although the origin for this anomaly is so far not established, many 
possible explanations for this effect have been considered in the literature, 
including modified gravity theories and models [3j, non-gravitational inter- 
actions H], extra-dimensional theories and models [5] and dark matter [6]. 

The main purpose of the present work is to investigate whether the Pi- 
oneer effect can be consistently obtained due to a background described by 
the ELA metric [ZllS], an ansatz that interpolates between the Schwarzschild 
(SC) metric [9] near massive bodies and the Robertson- Walker (RW) met- 
ric [TU] at spatial infinity, hence describing local matter in an expanding 
Universe |11] . The shift function of this metric can be fine-tuned through 
an exponential parameter allowing to describe intermediate scale gravita- 
tional effects [8l[12] which may be interpreted either as due to local cold 
dark matter (CDM) effects |13) or modified theories of gravity |14j . Al- 
though this fine-tuning is technically more challenging than semi-classical 
heuristic approaches, has the following advantages: 

• it is a covariant formulation allowing to deal with CDM within the 
framework of General Relativity; 

• in particular allows for direct covariant computations of gravitational 
accelerations and gravitational red-shifts (blue-shifts); 

• allows as well for QFT cross-section computations in curved back- 
grounds to include CDM effects (a metric is required to define the 
Dirac equation); 

• it allows for the description of intermediate scale CDM effects main- 
taining the local and global known physical properties of space-time. 

In the following we will employ two distinct coordinate systems, expand- 
ing coordinates, ri, for which the area of the sphere is ^ = Anrf such that 
the coordinate quantities correspond directly to measurable lengths and the 
standard Robertson- Walker coordinates, r, for which the area of the sphere 
is A = 47r a^r^ being suitable for direct computations of measurable red- 
shifts [H]. The map between both these coordinate systems is ri = ar and 
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the ansatz for the ELA metric is [71[8] 

= (1 - U){cdtf - a" (^^i^ + ^ (l - (1 - {cdt)j - a^dn' , 

(2) 

where a is the universe scale factor, H = a/a is the time dependent Hubble 
rate, U = 2GM/{c^ri) = 2GM/{c^ar) is the standard SC gravitational 
potential, M is the standard SC mass, c is the speed of light, G is the 
Newton gravitational constant and the solid angle line-elements are dO^ = 
rl{de'^ + sin^ ed^'^) and dJ^^ = r2(d^2 ^ Q^^2y 

This ansatz generalizes previous ansatz and solutions in the literature |15j 
allowing for local anisotropy [KJ and having the novelty of describing a 
space-time without singularities at the SC radius simultaneously maintain- 
ing as asymptotic limits the SC metric near the SC event horizon and the 
RW metric at spatial infinity. To ensure these properties of space-time the 
following limits must be obeyed [71[H] 

lim (1 — C/)° = : SC mass pole M is maintained , 

ri-->0 

lim a > 3 : space — time is singularity free at SC horizon , 

lim (1 — U)'^ = 1 : converges to RW metric . 

(3) 

In the following we consider as ansatz for the exponential parameter a a 
second order expansion in the gravitational field given by 

a(ri) = ao(ri) + a2(ri)(l - Uf , Qo(ri) = 3 , aiin ^ 0) < 0. (4) 

This ansatz allows for higher values of the corrections over intermediate 
scales than a first order expansion ansatz in the gravitational field maintain- 
ing the asymptotic limits expressed in (l3|). In particular, in the remaining of 
this work, we will heuristically fit the functional parameter 02 to the exper- 
imentally measured Pioneer effect which is fully accounted as a background 
blue-shift and compute a finite non-null mass for the extended correction to 
the background which may be interpreted as a CDM extended distribution. 

2 Describing the Pioneer Effect 

The Pioneer effect was measured using a two-way Doppler shift [T]. The 
standard gravitational red-shift as described by the SC metric is null, as for 
the standard cosmological red-shift is negligible as noted in being below 
the accuracy of the Pioneer measurements. For this reason the Pioneer effect 
is often attributed to a correction to the gravitational acceleration [3l[5l[6]. 
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Taking as central mass the Sun and neglecting the contribution of other 
bodies in the solar system, the corrections to the gravitational acceleration 
felt by a test mass travelling in the background given by the ELA metric ([2]) 
with respect to the General Relativity gravitational acceleration correspond- 
ing to a Ricci-flat background described by the SC metric is [8l[T2] 



K « ^ (l-C/0r(2(l-C/o)-(l + a)rp?7^ 



-2(1 + q)il- ^7©) 2-f + rp(l - Uq) ln(l - UQ)a') + 0{H^) 



^ +''P\ 

(5) 

where the primed quantities stand for differentiation with respect to the 
radial coordinate ri, rp stands for the geometrical (measurable) distance 
from the Sun to the Pioneer space-craft and Uq = 2GMq/ {c^Vp) stands 
for the standard SC gravitational potential of the Sun with Mq = 1.99 x 
10^0 Kg. 

In addition there is a correction to the gravitational red-shift. The mea- 
surable red-shift due to the background is obtained from the line element ([2]). 
For a light ray travelling approximately at a radial trajectory with respect 
to the Sun described by the four-vector ~ (a;/c, A;'', 0, 0) , from earth at 
Ti'^ = a^r''^^ ~ 1 AU to the Pioneer space-craft at r^'^ = oi r*-^-* = and 
back to earth at rj^^^ = a2 r^"^^ ~ 1 AU we obtain that the ratios between the 
radial wave numbers of the emitted radiation from earth A:|'q^ , the radiation 
received at the Pioneer space-craft k^^^^ and the radiation received back at 
Earth k''^^^ are 



,^ = ^(l-(l-f7)f-^) 



_ao l-Up-5 ^(2) _ ai l-Ue x - ^'^p 
fc^Q^ ai 1 -Ue ' A;^^^ 02 I - Up + 6 

(6) 

where Up and Ue are the gravitational potential of the Sun evaluated at the 
Pioneer space-craft and at Earth, and ao, ai and 02 stand for the Universe 
scale factor a = a{t) evaluated at the time of emission of the radiation 
from Earth, the time of receiving and re-emission at the Pioneer space-craft 
and the time of receiving it back at Earth, respectively. Further noting 
that expanding the Universe scale factor to first order in the Hubble rate we 
obtain that ai = ao{l+iH/c)rp)+OiH^) and 02 = aoil+2{H/c)rp)+0{H'^), 
the measurable red-shift (or blue-shift, if negative) due to the background 
for a two way Doppler measurement from Earth to Pioneer and back is 



Au _ ^("i) ^{2) 1 1 - Up 



1 + 0{H^) (7) 



Hence the contribution to the Pioneer non-physical acceleration due to the 
background red-shift is obtained by differentiating the standard two way 
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Doppler shift = —2fp/c 

where the prime denotes differentiation with respect to the radial geometric 
distance (the coordinate ri), Vp = rp for the radial geometric velocity of 
the space-craft (obtained from the differentiation dt = dtvi ) and we have 
neglected relativistic corrections (7 ~ 1 for fp <C c). 

Hence, if the Pioneer acceleration is totally due to the background en- 
coded in the metric ([2]), it is the sum of the contribution ^ and ([8]) 

ap = 5rp + 5r^°PP'- . (9) 

Assuming ao(^i) = 3 and an approximately linear functional exponent 
<^2{fi) = 0^2'^^ + a2^^ ri for radial distances above ri = 20 AU the Pioneer 
data is fitted, within the error bars, for the values of the constant coefficients 
= -1.23x10^ anda^^^ = -7.12 x 10"^ Fits within half of the error 
bars are also straight forwardly obtained as pictured in figure [TJ Within the 
range rp G [20, 70] AU the dominant contribution to ap is the background 
blue-shift, ^fDoppier ^ — 10~^^ ms"^. As for the value of the physical accel- 
eration, it is 4 orders of magnitude below this value, 5fp ~ +10"^'^ ms~^ 
(outwards from the Sun). 

(a) Pioneer acceleration witliiii tlie measure nic ill error bai'.'; (b) Pioneer acceleration within half of the measurement error bars 



20 30 40 50 60 70 20 30 40 50 60 70 

Figure 1: Fit of tKe^linear coefficient a2(^p) = <^2^^ + '^2^'* rp''io^ the Pioneer ac- 
celeration Op (QJj; (a) within the measurement error bars, = -1.225 x 10^ 
and a^'' — —7.118 x 10^"^ m^^; (b) within half of the measurement error bars, 
a^°^ = -9.459 X 10« and a^^^ = -7.067 x IQ-^ m-\ 



3 Compatibility with planetary motion 

Assuming the linear fit just discussed, within the Solar system, we obtain 
that the physical acceleration is positive above ri = 12.097 yl?7 increasing 
monotonically (at ri = 70 AU its value is Svp = 8.626 x 10"^*^ ms~^) and it is 
negative below ri = 12.097 ylf/ overcoming the Newton acceleration below 
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ri = 0.890 AU. These results are clearly not verified experimentally, the 
orbital motion within the solar system is driven and approximately described 
to a very good accm'acy by the Newton acceleration. This problem is straight 
forwardly solved by considering a cut-ofF at ri = 20 AU such that, below 
ri = 20 AU, this coefficient is approximately null. Hence we obtain the 
following heuristic ansatz 



with < a2 ^1 being an arbitrarily small non-null positive constant 
that ensures that the limits (131) are obeyed [7], 4^°^ = -1-225 x 10^ and 
02"^^ = —7.118 X 10^'^ m~^. We have also considered an upper cut-off Ruax 
that will be discussed in the next section. Although not a very elegant so- 
lution this heuristic ansatz allows to describe the Pioneer effect maintaining 
compatibility with the planetary motion within the solar system. 

Nevertheless the effects of this correction to the gravitational acceleration 
should be measurable at least for Pluto and Neptune: as an example, for 
the ansatz ()10p . if this correction is fully accounted by an orbital radius 
variation we would obtain, respectively, f orb/ f orb ~ 2.62 x 10~^^ 
and rorb/forb ~ 1-21 X 10^^^ which is well above the experimentally 
measured orbital radius variation for inner planets in the Solar system such 
as Venus and Mars, \rorb/'''orb\ < 10~^^ [TT]- However further analytical 
or numerical results are required to exactly predict the effects on planetary 
motion due to this correction [8] . 

4 Mass-Energy density and Possible Dark Matter 
Interpretation 

So far we have not analysed neither the mass-energy density neither the 
equation of state for the ELA metric ([2]). The mass-energy density (p) and 




(10) 
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anisotropic pressures {pr and pQ = p^p) are 



P 



1 - Uf-^ (a [/ + (1 - [/) (3 + n log (1 - U) dr,a)) , 

(2(1 + q) - (1 - U)"^ {aU 

+(l-[/)(3 + ri log{l-U)dr,a) 



Pr 



Pe 



8ttG 



levrG 



a-3 

[i-u)— 



a 



-2(1 - U)"^ ( 3 



4(l + g)(l-C/(l-- 



3[/ ( 2 - [/ ( 1 - I 



a 



1- 3 ) + n^7(i-f/)5., 



+ri(l - U) log(l - [/) ((1 + q)dr^a -2U{l-U)^ adr.a 
(l-U)"^ (6dr,a + n log(l - U){dr,af + nd^ 



a ^ / 

(11) 

We note that, for the hnear ansatz with one single branch for all values 
of the radial coordinate discussed in section EJ 02 = 02^ ''^ + 0^2^'^^ the 
mass-energy is negative below r\ = 8.06 AC/ hence violating causality. The 
branch ansatz (jlOp also solves this problem, for this ansatz p > for all 
values of the radial coordinate above the SC radius r\ > ri se = 2GM/cP. 

When the limits ([3]) are obeyed [7], at ri = ri.sC) P coincides with the 
mass-energy density for Ricci flat space-times, being null p(ri — ri se) = 0, 
and at spatial infinity it asymptotically converges to the RW mass-energy 
density p{ri — )• +00) = pRw = 3/i^/(87rG), as well as Pr{ri — )• +00) = 
Peii"! — ^ +00) = Pip{ri — )• +00) = prw = (2g— 1) /i^/ (Svr G) such that at spa- 
tial infinity the RW equation of state is obtained ujrw = (2g— 1)/3 ~ —0.726 
(this value is computed within the ACDM model assuming that the pressure 
of the background is given only by the cosmological constant [8]). Hence 
this mass-energy distribution has been interpreted in [7j, for a constant func- 
tional exponent a as a local anisotropic extended correction to the global 
isotropic cosmological background due to a local point mass vanishing at 
spatial infinity. In addition, when considering a radial coordinate depen- 
dence on the functional parameter 02, it has been put forward in [8|[T2] 
that the ELA metric ([5]) may also encode the corrections to the expanding 
background due to the existence of a local CDM distribution. If this is the 
case the total CDM mass can be computed corresponding to the mass of 
such extended corrections to the cosmological background. Specifically the 
CDM mass within a shell of radius Ri for a continuous ansatz is 

Mcdm{Ri) = 47r 1"^' rlip - pRw)dn = ^ ((1 " U{Ri)T^'^'^ - l) i?? • 

(12) 
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As long as the limits ([3]) are obeyed Mcdm{Ri — )■ 0) = 0. For convergence of 
this quantity at spatial infinity more severe conditions than the limits ([3]) are 
required, in particular it can be fine-tuned either to be null (e.g. for a{ri — )• 
+oo) = const/r", n > 2) or to a multiple of the central mass (e.g. for a(ri — )• 
+oo) = const/rf we obtain Mcdm(Ri +oo) = -H'^ const Mq/{2G)). How- 
ever for the ansatz (fTO|) for 02, Mcdm{Ri — ^ +00) = +00 and there is a 
discontinuity at r = 20 AU contributing a finite non-null contribution to the 
total mass. 

Further recalling that in average the ratio of CDM to baryonic matter 
is ficdm/^fe = 5 [E], if the EL A metric corrections are to be interpreted as 
due to a CDM distribution, each individual body should have a non-null 
finite value for the quantity Mcdm contributing to 0,cdm- This is straight 
forwardly achieved by considering a non-continuous upper cut-off^ for a at 
some value of the radial coordinate ri = i?Max such that above this value 
a = ao = 02 = and the metric for large radii coincides with the origi- 
nal isotropic McVittie metric [15] for which p = pRw- The CDM mass of 
this extended configuration is given by McDM{RMa.K) + -^cj:>a/(20~ AU) — 
McdmC^O'^ AU) McDMiRuax) and the cut-off i?Max corresponds to the 
edge of the CDM extended distribution. For the Sun, assuming the ratio 
Mcdm/Mq = 5, we obtain i?Max = 38317 AU. We remark that this upper 
cut-off also ensures that isotropy is exactly retrieved above -RMax- 

As for the equation of state for the ELA metric, we have already notice 
that at spatial infinity it coincides with the standard RW equation, hence 
neither the local baryonic matter correction to the background, neither the 
local CDM matter correction to the background modify the global equa- 
tion of state: statistically the respective pressures are null ph = Pcdm = 0, 
only the masses contribute to the global background densities fift and ricdm- 
However we do note that locally there are deviation from the global densi- 
ties, pressures and equation of state. For baryonic matter with a constant 
exponent a, the equation of state was analysed in detail in [7]. As for the 
extended distribution in the range r G [20 AU, -Rmax] corresponding to the 
ansatz ([10]), we obtain at r = 20AU, p = 2.71x10^ PRW,Pr = 3.74xl0^pw 
and Pe = P<fi = 2.74 x lO^p/jvK such that = Pr/p = — 1 + 6.81 x 10^® 
and u)0 = oj^p = pe/p = —7.33. These values change monotonically up 
to r = iimax where we obtain p = 1.36 x lO^p/jvK, pr = 1.87 x lO^p/jvy 
and Pe = Pip = 1-87 x l^'^ prw such that a;,. = — 1 + 7.44 x 10~^ and 
uq = oJip = —1 + 2.18 X 10~^. We further remark that the values of the equa- 
tion of state, —1 < Wr, We < —0.726, indicate that this construction may be 
compatible with a locally anisotropic configuration of scalar fields [T9] . 

Again we note that, although a clearly unelegant heuristic solution not 
derived from first principles, this construction allows for the interpretation 
of the Pioneer effect as due to dark matter, as well as, generally, allowing 
for a covariant description of local dark matter maintaining compatibility 
with today's global parameters of the universe. 
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5 Conclusions 



In this work we have considered a hnear ansatz for the functional parameter 
02 of the EL A metric ([2]) which allows to describe the Pioneer effect. The 
main contribution is due to the background gravitational blue-shift such that 
the physical acceleration contribution to the Pioneer acceleration is only 
of order ~ 10~^^ms~^, hence 4 orders of magnitude below the measured 
Pioneer acceleration 

Within this framework are also evaluated the orbital radius variation r jr 
for Pluto and Neptune. Orbital motion corrections to these planets are a 
common consequence of any gravity theory or model that allows to describe 
the Pioneer effect |20] . However we remark that the construction suggested 
in this work allows to significantly decrease the magnitude of such corrections 
due to the Pioneer effect being mostly attributed to a background blue-shift 
instead of a correction to the gravitational acceleration. 

In addition we also discuss the possibility of interpreting this effect as due 
to a CDM extended distribution from r\ = 20 AU up to a cut-off ri = i^Max 
with a computable fine-tuneable total mass Mqdm- 
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